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I . Abstract. We are reporting on a search for flavour-changing neutral current decays D" — > 

f^ l' M^/^~ &iid D'' — >■ e^e~, and for lepton-flavour violating decays D^ — >■ e fj,^ , the measurement 

Qj , of Df^ — )■ KgiT^ and Df^ -> KgK'^ branching fractions, and the search for CP violation in 

D+j -> K'^Tv+ and ^3+, -i> -R's-K'^ decays. The analyses are based on 600 fb"^ to 700 fl3~^ of 
data collected in e^e^ collisions at the centre-of-mass (CM) energy of the T(4S') resonance and 
60 MeV below by the Belle detector at the KEKB collider. 
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• I 1. Search for D" —?■ i^l decays 

f^ ■ The flavour-changing neutral current (FCNC) decays D^ — )• e~^e~ and D^ -^ fJ'^l-'-^ [1] are highly 

f^ I suppressed in the standard model (SM) by the Glashow-Iliopoulos-Maiani mechanism [2]. With 

the inclusion of long distance contributions the branching fractions can reach values of around 

10~^^ [3], these predictions are orders of magnitude below the current experimental sensitivity. 

The lepton-flavour violating (LFV) decays D^ -^ e^ji^ are forbidden in the SM. In certain new 

;^ ■ physics scenarios, FCNC branching fractions can be enhanced by many orders of magnitude and 

5^ , FLV decays may become possible [4]. For example, so far unobserved lepto quarks could enhance 

B (L>° -f /i+/x") to 8 X 10-'^ [5] . 

Using 660 fb^ of data we searched for the decays D^ — )• fi^^~ , D^ — )• e~^e~ and D^ — > e /x^. 
We use -D" mesons from the decays Z?*"*" — )• D^nf with a characteristic low momentum pion, 
since this considerably improves the purity of the reconstructed samples. We normalise the 
sensitivity of our search to topologically similar D^ — )• vr^vr^ decays; this cancels various 
systematic uncertainties. The signal efficiencies £« and EvrTr are evaluated using signal Monte 
Carlo simulation. 

In order to avoid biases, a blind analysis technique has been adopted. As the D^ -^ i^£~ 
decays are not expected to be observed at the current sensitivity, we maximise the figure-of- 
merit, T = e^/A^uL, where eu is the efficiency for detecting D^ — )• i~^£~ decays, and A'ul is 
the Poisson average of Feldman-Cousins 90% confidence level upper limits on the number of 
observed signal events that would be obtained with the expected background and no signal [6]. 

The background events can be grouped into two categories: (1) a smooth combinatorial 
background, and (2) a peaking background from the misidentification of D^ — >■ 7r''~7r~ decays. 
To estimate the number of combinatorial background events in the signal region, the sideband 
region is used. The peaking background in the signal region due to misidentification of 
D^ — > vr+vr" decays is estimated from the reconstructed D^ — > vr+vr^ decays found in data 



and the misidentification probability measured in data using D*'^ 
decays, binned in particle momentum p and cosine of polar angle. 
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Figure 1. The dilepton invariant mass distributions 
for a) D^ -^ fJ'^fJ'^, b) D^ -^ e+e^ and c) 



D -^ e ji^ . The dashed vertical lines indicate the 
optimised signal window. Superimposed on the data 
(open histograms) are the estimated distribution 
for combinatorial background (filled histogram), the 



misidentification of D 







vr ' TT 



(cross-hatched 



histogram) , and the signal if the branching fractions 
were equal to the 90% confidence level upper limit 
(single hatched histogram). 



The invariant mass distributions after applying the optimised event selection criteria are 
shown in Figure 1. In the signal region we find two candidates in the D^ — > fi'^fj,^, zero 
candidates in the D^ -^ e~^e~ and three candidates in the D^ -^ e /i^ decay mode; the 
yields are consistent with the estimated background of 3.1 ±0.1, 1.7=b0.2, and 2.6 it 0.2 events 
respectively. A binned maximum likelihood fit is used to determine the yield of D^ -^ it~^tt~ 
candidates for the normalisation. Finally, the branching fraction upper limits (UL) are calculated 
using the program pole.f, which extends the Feldman-Cousins method by the inclusion of 
systematic uncertainties [7]. The upper limits on the branching fractions at the 90% confidence 
level are found to be B [D^ -^ /U+M") < 1-4 x lO^^ B {D^ -^ e+e") < 7.9 x 10"^ and 
B (^D^ —7- e^/i^) < 2.6 X 10^^ [8]. Our results improve the current limits by a factor of 9 for 
-D" — )■ ijl^ iJi~ decay, by a factor of 15 for D^ — )■ e^e~ decay and by a factor of 3 for D^ -^ e^/i^ 
decay [9]. In 2008 the CDF collaboration reported a preliminary result on the UL for the 
D^ -^ l-i^l-i^ branching fraction [10]; our result is lower by a factor of 3 and strongly disfavours 
a leptoquark contribution [5] as the explanation for the anomaly in the measured D^ 
width [11]. 
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2. Branching fraction measurement of Df^ 
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K'^K+ decays 



Decays of charmed mesons play an important role in understanding the sources of SU(3) flavour 
symmetry breaking [12]. For D'^ decays, the branching ratio B{D'^ —^ K K~^)/B{D^ — > K vr"*") 
deviates from the naive tan^ 9c expectation [9], due to the destructive interference between 

colour-favoured and colour-suppressed amplitudes in D^ -^ K n^ [13]. However, converting 
experimental measurements of D decays that include Kg branching ratios to those involving 

K^ or K is not straightforward due to the interference between the doubly Cabibbo-suppressed 
(DCS) and Cabibbo-favoured (CF) decay modes where the interference phase is unknown [14, 15]. 
Based on a data sample of 605 fb~^ we measured the D~^ — ?> K'^K^ and D^ -^ K'^tt'^ 
branching ratios with respect to the corresponding Cabibbo-favoured modes. The invariant 
mass distributions of the selected events are shown in Figure 2. The results are B{D~^ —?■ 



KlK^ 



KlK+)/B{D+ -^ K'^TT^) = 0.1899±0.0011±0.0022 and B{Df -^ K^g7r+)/B{D+ 
= 0.0803ib0.0024ib 0.0019, where the first uncertainties are statistical and the second are 
systematic [16]. Using the world average values of CF decay rates [9], we obtain the branching 
fractions B{D+ -^ K^s^+) = (2.75±0.08)xl0"3 and B{D+ -^ KIti+) = (1.20±0.09)xl0-3 
where the uncertainties are the sum in quadrature of statistical and systematic errors. These 
are consistent with the present world averages [9] and are the most precise measurements up to 
now. The ratio B{D+ -^ K'^K+)/B{D+ -^ K'^tt^) = 2.29±0.18 may be due to SU(3) flavour 
breaking and/or different final-state interactions in D^ and Df decays. 
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Figure 2. Invariant mass distribution of selected KgK^ pairs LEFT and K^t:^ pairs RIGHT. 
Points with error bars show the data and histograms show the results of the fits. Signal, 
D|^N -^ Kgh'^ background {h = K ov tt), and random combinatorial background components 

are also shown. The inset is an enlarged view of the Df region. 



3. Search for CP violation in D, 



Kh' 



,. —r ±^gn and DT^x -^ KgK^ decays 
Another important aspect of such decays is the violation of the combined Charge-conjugation 
and Parity symmetries {CP). In the SM, the charmed particle processes for which a significant 
non-vanishing CP violation is expected are singly Cabibbo-suppressed (SCS) decays in which 
there is both interference between two different decay amplitudes and a strong phase shift from 
final state interactions. In the SM, CP violation in SCS charmed meson decays is predicted to 
occur at the level of 0(0. 1)% or lower [17]. 

Based on a data sample of 673 fb^^ we determine the CP violating asymmetry Aqp by 
measuring the signal yield asymmetry Arcc = (^rcc — -^rcc)/(-^rcc + -^rec) where N^(.c{N rec) is 
the number of reconstructed decays of DfJDV.). The measured asymmetry in this equation 
includes two contributions other than Acp- One is the forward-backward asymmetry (Aps) due 
to 7* — Z^ interference in e~^e~ — )■ cc and the other is a detection efficiency asymmetry between 
positively and negatively charged tracks A^ = (e^ — e^)/(e^ + £ )) where e^(e^) is the efficiency 
for K~^{K^) or 7r"'"(7r~) meson and h denotes K or n. Since K^ mesons are reconstructed from 
a 7r"'"7r~ pair, there is no detection asymmetry other than A^. The signal yield asymmetry can 
therefore be expressed as Arec = Acp + App + A^. 

To correct for the asymmetries other than Ac p., we use reconstructed samples of Df — )■ 07r+ 
and D^ -^ K~tt~^ decays and assume that Acp in CF decays is negligibly small compared to 



the current experimental sensitivity and that ApB is the same for all charmed mesons. We 



reconstruct 6 mesons via their 



K^K decays. The measured asymmetry for D^ 



is the sum of ApB and A^. Hence one can extract the Acp value for the if^7r+ final state by 
subtracting the measured asymmetry for Df 



ivr"'" from that for D 



is) 



i^07r^ 



The method for the measurement of Acp in the KgK^ final states is different from that for 
the KgiT^ final states. The App and A'^ components in A^ec are directly obtained from the 
Df — )■ (j)TT~^ sample, but there is no corresponding large statistics decay mode that can be used 
to directly measure the App and A^ components in ^rec- Thus, to correct the reconstructed 
asymmetry in the KgK~^ final states, we use samples of D^ — )■ K~tt~^ as well as D'^ -^ cpir^ 
decays. The value ^rec — A^ includes not only an Acp component but also an App component. 
Since Acp is independent of all kinematic variables, while App is an odd function of cos6^^^, 

we can deduce both by addition/subtraction in bins of cos 9. Figure 3 shows the results. 
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Figure 3. Measured Acp and App values for 
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as a function of | cos 6^^^ \ . The 
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dashed curves show the leading-order prediction 
for Af^. 



No evidence for CP violation has been observed [18]. Our results are Af-,p ^ = (—0.71 it 

0.19 ± 0.20)%, ^J^J,^-^°^^ = (+5.45 ± 2.50 ± 0.33)%, A^p~^^°^^^ = (-0.16 ± 0.58 ± 0.25)%, 

and Aqp "^ = (+0.12 ± 0.36 ± 0.22)%. They are consistent with the SM predictions and 
provide the most stringent constraints up to now on models beyond the SM [14]. 
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